Statistics of work performed on a forced quantum oscillator 
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Various aspects of the statistics of work performed by an external classical force on a quantum 
mechanical system are elucidated for a driven harmonic oscillator. In this special case two parameters 
are introduced that are sufficient to completely characterize the force protocol. Explicit results for 
the characteristic function of work and the respective probability distribution are provided and 
discussed for three different types of initial states of the oscillator: microcanonical, canonical and 
coherent states. Depending on the choice of the initial state the probability distributions of the 
performed work may grossly differ. This result in particular holds also true for identical force 
protocols. General fluctuation and work theorems holding for microcanonical and canonical initial 
states are confirmed. 
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I. INTRODUCTION 

During the last decade various fluctuation and work 
theorems [1, 2] have been formulated and discussed. 
They inter alia characterize the full nonlinear response 
of a system under the action of a time dependent force 
[3, 4]. These theorems have been derived and experi- 
mentally confirmed primarily for classical systems [5-7]. 
Quantum mechanical generalizations were proposed re- 
cently [8-15]. 

Conceptual problems though arise in the context of 
quantum mechanics if one tries to generalize those clas- 
sical relations that require for example the specification 
of a system's trajectory extending over some interval of 
time, or the simultaneous measurement of noncommut- 
ing observables. For example, the measurement of work 
performed by an external force on an otherwise isolated 
system may be accomplished in the framework of classi- 
cal physics in principle in two different ways. The first 
method is based on two measurements of the energy, one 
at the beginning and the second at the end of the consid- 
ered process. This method becomes unreliable in prac- 
tice if the system is large and the work performed on the 
system is negligibly small compared to the total energy 
of the system. Such a situation typically arises if the 
system of interest, on which the force exclusively acts, 
interacts with its environment. In order to retain an iso- 
lated system, the large system made of the open system 
and its environment must be considered. Again, the work 
performed on the system results as the difference of the 
energies of the total system, which both may be very 
large. 

For classical systems, this unfortunate situation can 
be circumvented by a second method, by monitoring the 
state of the relevant small system during the time when 
the force is acting. Having this information at hand one 
can determine the work by integrating the power sup- 
plied to the system at each instant of time. The respec- 
tive power can be inferred from the registered state of 
the system and the known force protocol. In a quan- 
tum system a continuous measurement of even a single 



observable would strongly influence and possibly mani- 
festly distort the system's dynamics. Apparently, only 
the first method of two energy measurements is feasible, 
at least in principle, in the quantum context. 

An alternative method based on a continuous mon- 
itoring has recently been suggested by Esposito and 
Mukamel [11] for open quantum systems described by 
Markovian quantum master equations. There the dy- 
namics of the density matrix is mapped onto a classical 
rate process for which known fluctuation theorems can 
be applied [16]. This provides an interesting formal ap- 
proach but its physical meaning has remained unclear 
[11]. Moreover, this approach is restricted to open sys- 
tems that only weakly interact with their respective en- 
vironments. 

In the present paper the distribution of work is dis- 
cussed for the exactly solvable system of a driven har- 
monic oscillator [15, 17]. In this case, the distribution of 
work is discrete. We provide formal expressions for this 
distribution and its corresponding characteristic function 
which are valid for all initial states of the system as well as 
for all possible kinds of force protocols. In particular, we 
determine the characteristic functions and distributions 
of the work for microcanonical, canonical and coherent 
initial states which lead to qualitatively different work 
distributions. 

The paper is organized as follows. In Sect. II we re- 
view the general form of the characteristic function of 
work performed on a system in terms of a correlation 
function of the exponentiated Hamiltonians at the initial 
and final time of the force protocol. We prove that this 
particular expression indeed always represents a charac- 
teristic function, i.e. the Fourier transform of a proba- 
bility density. Sect. Ill presents various fluctuation and 
work theorems for canonical and microcanonical initial 
states. In Sect. IV general expressions for the character- 
istic function and the corresponding probability distribu- 
tion of work are derived for a driven harmonic oscillator. 
Moreover, the expressions for the first four cumulants are 
derived. The dependence of the work distribution on the 
force protocol for microcanonical, canonical and coher- 
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ent initial states as well as its dependence on the specific 
parameters of these initial states are investigated. 



II. CHARACTERISTIC FUNCTION OF WORK 

The response of a quantum system on a perturbation 
by a classical, external force can be characterized by the 
change of energies contained in the total system. The 
energy as an observable coincides with its Hamiltonian 
H(t) of the total system. It includes the external force 
and therefore depends on time. We will consider the dy- 
namics of the system only within a finite window of time 
[to,tf] during which the force is acting in a prescribed 
way, resulting in a protocol of Hamiltonians which is de- 
noted by {H{t)} tfto . Apart from the action of the exter- 
nal force the system is assumed to be closed. Its dynam- 
ics is consequently governed by a unitary time evolution 
Ut^oi which is the solution of the Schodinger equation 

ihdu t , t0 /dt = H(t)u t , t0 , 

u(t,t ) = i. u 

As explained in the introduction, the work w is measured 
as the difference of the energies of the system at the fi- 
nal and initial times tf and to- I n a single measurement 
the work is given by the difference of two eigenvalues 
e n (tf) and e m (io) of the Hamiltonians H(t) at the re- 
spective times tf and to, i.e. by w = e n (tf) — e m (to)- 
The inherent randomness of the outcome of a quantum 
measurement in general leads to a measured work that is 
random. A complete description of the statistical prop- 
erties of the work performed on the system is provided by 
the characteristic function G t0it/ (w) which presents the 
Fourier transform of the probability density of the work 
Pt t ,t (w), i.e. 

Gt ,t f (u)= J dwe iuw p tftt0 (w). (2) 

It can be expressed as a quantum correlation function 
of the exponentiated Hamiltonian at the initial and the 
final time [14], i.e. 

G t0 , tf (u) = (e iuH ^e- iuH ^} 

= Tre luHH ^^e- luH ^p{t Q ), 

where 

H H (t f ) = U+, to H(t f )U tf , to (4) 

denotes the Hamiltonian in the Heisenberg picture. The 
density matrix p{to) from the initial density matrix p(to) 
as a result of the measurement of the Hamiltonian H(to). 
It is given by 

P(t ) = Yl P n(to)p(to)Pn(t ), (5) 
n 



where the operators P n (to) denote the eigenprojection 
operators of the Hamiltonian at time to , which present a 
partition of the unity 

J2Pn(t ) = l. (6) 

k 

Before we apply the general expression (3) to a particular 
system and investigate its dependence on the initial state 
p(ta), we discuss three general properties of the correla- 
tion expression (3) which guarantee that the resulting 
function G t/jto (w) indeed always presents a proper char- 
acteristic function of a classical random variable w. This 
is the consequence of the three following properties: 

(i) G tft t a {u) is a continuous function of u. 

(ii) G t/ , to (u) is a positive definite function of u, i.e. for all 
integer numbers n, all real sequences u\, U2, • • • u n , and 
all complex numbers Zi, i = 1, 2 . . . n 

n 

Y, G t f ,t {ui-Ui')z*Zi' > (7) 

holds. Here, the asterisk z* denotes the complex conju- 
gate of Zi. 

(m) tVo(o) = i 

According to a theorem by Bochncr [18] the properties 
(i-iii) are necessary and sufficient conditions in order that 
the function G tf ,t { u ) is the Fourier transform of the 
probability measure of a random variable. In short, the 
first condition insures that, strictly speaking, the func- 
tion G tf ,t (u) is the Fourier transform of a measure, the 
second condition assures that this measure is positive and 
the third condition that it is normalized. Hence the cor- 
relation expression eq. (3) always defines a proper char- 
acteristic function. For a proof of the properties (i-iii) we 
refer the reader to the appendix A. 



III. CANONICAL AND MICROCANONICAL 
INITIAL STATES 

In experiments an external force is often applied on a 
system that initially is found in a thermodynamic equi- 
librium state. Depending on whether the system was in 
weak contact with a heat bath or was totally isolated 
from its environment, the initial state of the system is 
described either by a canonical or a microcanonical den- 
sity matrix. For both situations fluctuation and work 
theorems are known. We will shortly review these rela- 
tions. 



A. Work and fluctuation theorems for canonical 
initial states 

If the initial density matrix is canonical, i.e. if 

p(t o ) = Z- 1 {to)exp{-0H{to)}, (8) 
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where 



Z(t ) = Trcxp{-/?# (t )} - 



(9) 



denotes the partition function and F(to) the free energy, 
then [H(to), p(to)] — and the first measurement of the 
energy leaves the density matrix unchanged, such that 
p(to) = p(to). With cq. (3) this leads to the character- 
istic function of work for a canonical initial state which 
was derived in Ref. [12]. In this case, Gt f ,t (u) can be 
continued to an analytic function of u for all < ^su < [3 
[13]. For the particular value u = i(3 the characteristic 
function yields the mean value of the exponentiated work, 
(exp{— f3w}) and the correlation function expression (3) 
simplifies to the ratio of the partition functions at the 
times tf and to, resulting in the Jarzynski work theorem 

{e-O") = Z(t f )/Z(t ) = exp{-(3(F(t f ) F(t ))} , 

(10) 

where Z(t f ) = trexp{ -/?#(*/)} = exp{-0F(t f )}. 

Within the domain of analyticity S — {u\0 < 3u < 
/3} the characteristic functions for the original and the 
time reversed protocol are related to each other by the 
following formula, cf. [13] 

Gt f ,t (u) = ^G t0 , tf (-u + if3), (11) 

where Gt Q .t f {u) refers to processes under the time re- 
versed protocol {H{t)}t ,t f starting from the canonical 
state Z- 1 (t f )exp{-/3H(t f )}. An inverse Fourier trans- 
form leads to the Tasaki-Crooks fluctuation theorem, 
which relates the probability densities of work Pt f ,t ( w ) 
for a given protocol to the density of the work Pt ,tf{w) 
for the time reversed protocol. This theorem explicitly 
reads [13] 



PtfjtoH = Zjtf) j iw = e -f3(F(t f )-F 

Pto,t,{-w) Z{t ) 



(t )-w) 



(12) 



B. Fluctuation theorems for microcanonical initial 
states 

A system in a microcanonical state is described by the 
density matrix 



p{t )=w E l {t )8{H{t )-E), 



(13) 



where 



™ E (t ) = Tr 5(H(t ) -E)= cxp {S(E, t Q )/k B } (14) 

denotes the density of states as a function of the energy 
E of the system. The density of states can be expressed 
in terms of the entropy of the system Se(£o) provided 
the spectrum of the system Hamiltonian is sufficiently 
dense such that the density of states becomes a smooth 
function on a coarsened energy scale. The microcanonical 



density matrix commutes with the Hamiltonian H(to). 
Consequently, p(to) and p(to) coincide. 

The microcanonical quantum Crooks theorem assumes 
the form [14] 



Pt f ,t (E,w) 
Pt ,t f (E + w,-w) 



ZU E {to) 

S(E- 



exp 



j,t f )-S(E,t ) 



k B 



Y 

(15) 



Analogous to the canonical case it relates the probability 
density p to ,t {F 7 w) of work w, for a system starting in 
a microcanonical state with energy E, to the respective 
quantity for the time reversed process starting at energy 
E + w. This quantum theorem is formally identical to 
the respective classical theorem [19]. 

From the microcanonical Crooks theorem the proba- 
bility density relating to the time reversed process can 
be eliminated to yield the so-called entropy-from-work 
theorem [14], reading: 



^E f (tf) — J dw zu Ef - 



>{to)Pt f ,t { E f ~ w,w) . (16) 



This theorem allows one to determine the unknown den- 
sity of states of a system with Hamiltonian H{tf) from 
the known density of states of a reference system H(to) by 
means of the statistics of the work that is performed on 
the system in a process that leads from the reference sys- 
tem to the final system with unknown density of states. 
In the case of systems with a sufficiently smooth density 
of states the respective entropy can be determined. For 
further details see in Ref. [14]. 



IV. DRIVEN HARMONIC OSCILLATOR 

To illustrate these concepts we consider an example 
which allows the analytical construction of the probabil- 
ity of work. Specifically we consider a harmonic oscilla- 
tor on which a time dependent force acts during a finite 
interval of time. Its time evolution is governed by the 
Hamiltonian 



H(t) = hua+a + f*(t)a + f(t)a + , 



(17) 



where u denotes the angular frequency, and a + and a 
creation and annihilation operators, respectively, which 
obey the usual commutation relation, i.e. [a, a + ] = 1. 
The complex driving force f(t) allows for a coupling to 
position and/or momentum of the oscillator. We assume 
that f(t) vanishes for times t < to = 0. It is our aim 
to study the influence of the initial state p(to) on the 
statistics of work performed on the oscillator. The mea- 
surement of H(t ) = huja + a at time t = then yields 
the result Hujii with probability 



Pn = ( n |p(*o)|«)- 



(18) 



4 



Accordingly, the oscillator is found in the state 

p(to) = ^>r»<n| (19) 

n 

immediately after this measurement. Putting this den- 
sity matrix in the general expression for the characteristic 
function, eq. (3) one obtains 

G t/lt0 («) = j2p" e ~ iuhwn ( n \ eiuHH{t,) \ n ) ■ ( 2 °) 

n 

i 



For the driven harmonic oscillator the diagonal matrix 
element of the exponentiated Hamiltonian Hn(tf) can 
be determined [17]. For details see the Appendix B. 
With the expression (B13) for the matrix clement 
(n\ exp {iHh (tf)} \n) we find 



Gtf ,M =e »l/(*/)| 2 /(M exp (( e »^ _ 1} |z|2} J2Y,Pn 

= e H/(*/)| 2 /(M 



n=0 k=0 



|2(n-fc) 



k) (n — k)\ 

exp{(e^ - 1) \z\ 2 } f^PnU^sin 2 

n=0 ' ' 



g— iuhu>(n— k) f^iuhu _ j^j 2 ^ n ^ 



(21) 



where \f(tf)\ 2 /(hw) denotes a uniform shift of the spec- 
trum of the harmonic oscillator due to the presence of 
the external force, cf. eq. (B8), and 

1 [ t; 

z = — / dsf(s) exp{iujs} (22) 
hw J 

is a dimensionless functional of the driving force f(t), 
cf. eq. (B6). This dimensionless quantity vanishes in 
particular for all-quasi static forcings, i.e. if the force 
changes only very slowly with f(t) = g{t/tf) for tf — ► oo, 
where g(r) is a continuously differentiable function for 
t e [0, 1]. We hence call z(t) the rapidity parameter of 
the force protocol. Finally, L n (x) — J2k=o (1) ( — x) k /k\ 
denotes the Laguerre polynomial of order n [21]. 

Introducing the cumulant generating function K(u) = 
In G(u) one obtains the cumulants of work k n as the nth 
derivatives of K(u) with respect to u taken at u = [22], 
i.e. k n = (-i) n d n K(0)/du n . The first four cumulants 
become: 



ki = (w) 
_ \f(tf)\ 2 



fku 

k 2 = (w 2 ) - (w) 2 



+ hjj\z\ 



2{h^y\z\ 2 (a + a} 



(23) 
(24) 
(25) 



2) ' 

k 3 = (w 3 ) -3(w 2 )(w) + 2(w) 3 

= (M 3 M 2 . 

k 4 = (w 4 ) -4(w 3 )(w) -3(w 2 ) 2 
+ 12(w 2 )(w) 2 -6(w} 4 
= {Hujf \z\ 2 {1 + 4(a+a) + 6 [(a + a(a + a - 1)) 
-2(a+a) Q ]\z\ 2 }. (26) 



The odd cumulants of the work are independent of the 
initial preparation. The even cumulants depend on the 
factorial moments of the number operator a + a with re- 
spect to the initial state p(to) such as (a + a) = J2 n np n 
and {a + a(a + a — 1)) = J2 n n ( n ~ l)Pn> where p n is de- 
fined in eq. (18). Moreover, all cumulants apart from 
the first one vanish for forcings with z = 0. This holds 
true in particular for all quasi-static force characteristics. 
The underlying work probability density then shrinks to 
a delta function at w = \f(tf)\ 2 /(hu). 

In general, the work probability density follows 
from the characteristic function by means of an in- 
verse Fourier transformation. Rather than the char- 
acteristic function itself we first consider the function 
Q(u) = exp{-iu\f(tf)\ 2 /(huj)} G tf , to (u). Upon expand- 
ing exp {|z| 2 exp{iw?iaj}} into a series of powers of \z\ 2 
we obtain for Q(u) a Laurent series in the variable 
exp {iuhuj} . The inverse Fourier transformation is given 
by a series of delta functions 5(w — fuvr), with r e Z, 
with weights 

oc n 2k 

*=e- W ' E EEt- 1 )"^ 

m,n=0k=0 1=0 
| z |2(k+m) / n X / 2fc X 

x ^kT{k [i 6l+nhk+r , , 
\ / \ / , 27 

oo n min(fc+r,2fe) v ' 

= e " w2 EE E (-!) 2fc -^ 

n=0 k=0 1=0 
U|2(2fc+r-I) /„\/ 2 fc> 



{k + r-l)\k\ \kj V I 



The factor exp {— iu\f(tf)\ 2 /(tku)}, by which Q(u) has to 
be multiplied to yield Gt / ,t (u), gives rise to a constant 
shift such that the probability density of work performed 
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on a harmonic oscillator assumes the result 

\f(tf)\ 2 



Pt f ,o(w) — 2jgy S( w - (hujr + 



hut 



(28) 



In the next Section we will investigate the influence of 
the initial state on the statistics of the work. 



A. Distributions of work for different initial states 

As particular examples of initial states we will discuss 
microcanonical, canonical and coherent states. 



1. Microcanonical initial state 

For a microcanonical initial state with energy Hlotiq the 
density matrix becomes 



p(*o) = P(*o) = K)(«o|- 
The characteristic function then reads 
G™ o (n , u) = e iu l/(t/)l 2 /(M exp | ( e *ufi« 

l l l2 • 2 F UjJU 

x L no \ 4|zrsm^ 



(29) 



-1)1*1*} 



(30) 



and, accordingly, the probability q" lc (no) to find a change 
of energy by w = hwr + \f(tf)\ 2 /(hio) emerges as 



"0 



i(fc+r,2fc) 



V\ /2fc 



(-1) 



2fc-Z 



(fc + r-Z)!fe! 



(31) 



: '/Jl/ ' 



2(2fc+r-i) 



As expected from the behavior of the moments, all proba- 
bilities g™ c (no) with r ^ vanish for quasi-static forcing, 
i.e. if z — > 0. The dependence of q™ c (no) for n = and 
3 as well as for the eight lowest values of r on the pa- 
rameter z is displayed in Fig. 1. With increasing values 
of the rapidity parameter z the distribution is becoming 
broader. For the fixed value of \z\ = 2 the distribu- 
tion g™ c (no) is compared for the three initial states with 
no = , 10 and 30 in Fig. 2. With increasing value 
of no the distributions become broader. They develop a 
slightly asymmetric shape with higher peaks at negative 
values of r compared to those at positive r values. Be- 
tween these dominant peaks the probability still displays 
pronounced variations. 

For a harmonic oscillator, the microcanonical Crooks the- 
orem reduces to the relation q™ c (n) — gf^(n + r). One 
can show that this symmetry is fulfilled by the probabil- 
ities q™ c (n) given by eq. (31). 




FIG. 1: The probabilities q™ c (no) for two microcanonical 
initial states with no = (a) and no = 3 (b) are depicted 
for r = —3 ... 4, as functions of the rapidity parameter z in 
eq. (B6). In both cases the distribution collapses at r = for 
adiabatic forcing corresponding to \z\ =0 and broadens with 
increasing \z\. Obviously, when starting in the ground state 
the oscillator cannot deliver work whence the probability for 
negative r strictly vanishes. "Stimulated emission" becomes 
possible from an excited state at finite driving rapidity z lead- 
ing to nonzero probabilities q™ c (no) at negative values of r in 
panel (b). 



0.15 




0.05 



FIG. 2: (Color online) The probabilities g™ c (no) for a mi- 
crocanonical initial state with no = (circles) and no = 10 
(diamonds) and no = 30 (crosses) are compared for a fixed 
rapidity parameter \z\ = 2 and r = —22 ... 30. The lines serve 
as a guide for the eye. 
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2. Canonical initial state 



For a canonical density matrix 

p(to) = (1 - e -/3^) e -/ 3fiwa+a 



(32) 



the initial states are distributed according to p n 



-ft hum 



/(i- 



). This allows one to perform the sum 



over n in the expression for the characteristic function 
(21) in terms of the generating function of the Laguerre 
polynomials, cf. [21] yielding the expression 



G c t , ta {P,u) =ex P 



mil 

HlU 



M 2 -4|z| 2 Sin j^ /2) 

11 11 pprm i 



(33) 



Putting u — if} one finds that the two terms in the ex- 
ponent which are proportional to \z\ 2 cancel each other, 
such that one obtains 

(e Pw ) = G c tf>to (f},if3) = exp{-/?|/(M 2 /(M} ■ (34) 

The free energy difference of two oscillators with Hamil- 
tonians i?(to) = ftwa + a and H(tf) — hwa + a + f*(tf)a + 
f(t)a + ) each one staying in a canonical state at the 
temperature fj is given by AF = F(tf) — F(t ) = 
\f(tf)\ 2 /(hu>) in accordance with Jarzynski's work the- 
orem. 

The probability q^{f}) to find the work w = Hu>r + 
\f(tf)\ 2 /(ftio) if the system starts in a canonical state 
becomes 



(!-.-») 



oo n min(fc+r,2k) 

EE E (~^-" n 

n=0 k=Q 1=0 



\2(2k+r-l) 



(k 



iy.k 



(35) 



where f} = fttuv denotes the inverse dimensionless tem- 
perature of the initial state. The expression for <?£(/?) can 
be further simplified to read 



-|z| 2 coth(/3/2) e /3r/2 / 



sinh/3/2 



(36) 



where I v (x) denotes the modified Bessel function of first 
kind of order v [21]. For details of the derivation see 
Appendix C. 

Note that the following detailed balance like symmetry 
relation exists, 



-$r 



(37) 



relating the occurence of positive and negative work. In 
Fig. 3 the z dependence of <j£(/3) for f} = ln(4/3) is com- 
pared for a few small values of r. One finds that due 
to the average over the canonical initial distribution the 
multipeaked structure of the microcanonical distribution 
as a function of the rapidity parameter \z\ disappears 
and only a single peak remains for each value of r. The 



Qr (ln(4/3)) 




FIG. 3: The probabilities q£(/3 = ln(4/3)) for a canonical ini- 
tial state are displayed for r = —3 ... 4 as a function of the 
parameter z. For the sake of comparability the dimensionless 
inverse temperature is chosen such that the average energy in 
the initial state coincides with the energy 3hui of the micro- 
canonical state in Fig. 1 (b). 
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FIG. 4: (Color online) The probabilities <jf (/3) for a canonical 
initial state are displayed as functions of r for \z\ =2 and 
different values of the dimensionless inverse temperature f) = 
0.5 (boxes), 1 (circles) and 5 (crosses). The lines serve as a 
guide for the eye. 



temperature dependence of the work distribution is il- 
lustrated in Fig. 4. Finally, we verify the validity 
of the Tasaki-Crooks theorem (12) for a driven oscilla- 
tor. For this purpose we consider the probability density 
Pt .t f {—w) for the time reversed protocol. Since the ab- 
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solute values of the rapidity parameters coincide for the 
original and the time reversed protocol the probability 
density of work for the reversed protocols becomes 

Pt ,t f (-w)= £ q r 8(-w-(t^r-\lM^)\ , 

r— — oo v 



(38) 

where we took into account the overall shift of the spec- 
trum by the reversed protocol as — \f(tf)\/(ftw). Multi- 
plying both sides of eq. (38) with exp{— /3(AF — w)} = 
exp {— (3(\f(tf)\ 2 /(fkd) — w)} one obtains 

e -0 { AF- w)pto t/ { _ w) = J- e -/3(l/(*/)| 2 /(^)-) ? c (/ 9) 



x 6 I -w - (fuvr - 



r 

= Pt f ,t (w), 



(39) 



in accordance with the Tasaki-Crooks theorem (12). 



3. Coherent initial state 

An oscillator prepared in a coherent state |a) is de- 
scribed by the density matrix 



p(t ) = \a)(a\ 



(40) 



where 



\a) 



aa + -\-aa 



(41) 



and |0) is the normalized ground state of the oscilla- 
tor satisfying a|0) = 0. Note that the coherent state 
density matrix does not commute with the Hamilto- 
nian H(to). The measurement of H(to) modifies the co- 
herent state (40) by projecting it onto the eigenstates 
\n) = (a+)"/\/n!|0) of this Hamiltonian leading to 



p(to) 



|2n 



\n)(n\ 



(42) 



This implies a Poissonian distribution of the respective 
energy eigenvalues fkon 



a 



In 



III 



-I«| 2 



(43) 



which yields for the characteristic function of work (21) 
a closed expression of the form 



+ z 



(e ih0JU - 1) | 



x Jo 4 



huju 



az sin ■ 



(44) 



where Jq(x) is the Bessel function of order zero, cf. 
Ref. [21]. For the probability qr S (a) of work one obtains 
with eq. (27) 



-N 



OC 



\z\ 2m \az\ 2 \ m - r 
rn\(\m-r\\f 



-1F2 



\m — r\ + — ; |m — r\ + l,2\m— r\ + 1; — 4\az\ 



(45) 



where iF 2 (a; b, x; x) denotes a generalized hypergeomet- 
ric function [21]. For details of the derivation see the 
Appendix D. 

The dependence of the probabilities q^ s (a) on the ra- 
pidity parameter |z| is illustrated in Fig. 5 for r values 
ranging from —10 to 20. Increasing values of z lead to a 
broadening of the distribution and also to a shift towards 
larger values of r, see also panel (b) of Fig. 6. This is in 
accordance with eq. (23) and (24) for the first two cumu- 
lants of the work, which both increase with |z| 2 . Panel 
(a) of Fig. 6 shows the dependence of the probabilities 
qr S (a) on the parameters a. Increasing a also leads to a 
broadening of the work distribution without influencing 
its mean value, cf. also eq. (23). 



In Fig. 7 the probabilities q r are depicted for differ- 
ent initial states. In one case the oscillator is initially 
prepared in a canonical state at inverse dimensionless 
temperature (3 = (3fno = 0.1. In the other case, the 
oscillator stays in a coherent state \a), where the ab- 
solute value of I a I is chosen such that the mean exci- 
tation number is the same for both states, i.e. |a| 2 = 
exp{— Phuj}/(1 — exp{— (3hoj}). For f3hu) — 0.1 one finds 
I a 1 2 w 9.51. The two oscillators then are subjected to 
protocols with the same rapidity parameter |z| = 2. Ac- 
cording to eqs. (23) and (24) the first two moments of 
the work performed on the oscillator coincide. Yet the 
distribution of weight factors q^(/3) and q^. s (a) distinctly 
differ. Whereas the distribution is pronouncedly bimodal 
in case of the coherent state, it is unimodal for the canon- 
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FIG. 5: The probabilities qr S (ct) for a coherent state with 
parameter |q| = 3 are displayed for r = —10 ... 20 as functions 
of z. 



0.2 



0.15 



0.1 



0.05 



■ (a) 




0.1 ■ 






1 






10 




— 1 ™H»»»»»»fr< 





-10 



10 



20 



0.25 



0.15 



o< 0.1 



0.05 - 



1 1 1 1 

■ (b) . 


1 1 1 


i ■ 

= i ■ 








= 2 " 








= 4 ■ 











10 



20 



30 



FIG. 6: (Color online) The distribution of work performed 
on an oscillator which initially is prepared in a coherent state 
| a) for different values of a in panel (a) and of the rapidity 
parameter z in panel (b). In panel (a) the rapidity param- 
eter has the value \z\ = 2. In panel (b) the coherent state 
parameter has the value \a\ 2 = 1. 



ical state. The weight factors q^{(3) almost perfectly fall 
onto the Gaussian probability density which has the same 
first two moments as the discrete distribution given by 
q r . 



FIG. 7: (color online) The distribution of work is compared 
for a canonical and a coherent initial state subject to the same 
force protocol with rapidity parameter z — 2. With f3hui = 0.1 
and \a\ 2 ~ 9.51 the expectation values of the energies agree 
in the two initial states such that according to eqs. (23) and 
(24) the first and second moments of the work also coincide. 
Still the distributions of work grossly differ from each other. 



V. CONCLUSIONS 

In this work we studied the statistics of work per- 
formed on an externally driven quantum mechanical os- 
cillator by means of a correlation function expression 
of the work. We demonstrated that this particular ex- 
pression indeed always represents a proper characteristic 
function of a random variable, which is the performed 
work in the present context. The proof given here is 
based on Bochner's theorem. It holds for general quan- 
tum mechanical systems, not only for harmonic oscilla- 
tors. 

The considered force linearly couples to the position 
and momentum of the oscillator. It may describe the 
influence of an electric field on charged particles in a 
parabolic trap or the external forcing of a single elec- 
tromagnetic cavity mode. For this type of additive forc- 
ing, the frequency of the oscillator remains unchanged 
and therefore the level spacing of the eigenvalues of the 
Hamiltonian is not influenced by the force. The spectrum 
is only shifted as a whole. As a consequence the work 
performed on the oscillator is, as a positive or negative 
integer of the level spacing, a discrete random variable. 
We determined the first few cumulents of the work for 
arbitrary force protocols and initial states. A comple- 
mentary study for a parametrically forced oscillator was 
recently performed by Deffner and Lutz [15]. 

It turns out that for the harmonic oscillator the statis- 
tics of work depends on the force protocol {f(t)} tf .t 
only through two real parameters, which are (i) the 
shift of the spectrum, given by L(tf) = \f(tf)\ 2 /(huj), 
and (ii) the absolute value of the dimensionless quantity 
z = J t f f(s)exp{iujs}. This parameter vanishes for all 
quasi-static processes and therefore presents a measure 
of the rapidity of the force protocol. While the presence 
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of L(tf) only causes an overall shift of the possible values 
of the work, the rapidity parameter \z\ also influences its 
distribution. Typically, the distributions move towards 
larger values of work w and become broader with increas- 
ing rapidity \z\. 

We also demonstrated that different initial states of 
the system such as microcanonical, canonical or coher- 
ent states, have a large influence on the work statistics. 
We further note that two different initial density ma- 
trices with the same diagonal elements with respect to 
the energy eigenbasis of the Hamiltonian H(t ) lead to 
identical work distributions even though the two density 
matrices may be very different in other respects. For ex- 
ample, the coherent pure state |a;)(a| and the mixed state 
exp{ — H 2 "} J2 n M 2 / n 'l n )( n l cannot be distinguished by 
means of their respective work statistics. This statistics 
is also insensitive to the phase of a coherent state. 
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APPENDIX A: PROOF OF THE PROPERTIES 

OF G t/ , t0 (u,v) 

We prove that the conditions of Bochner's theorem are 
fulfilled, and consequently G tf ,t { u ) is a proper charac- 
teristic function. 

Proof of property (i): G t , to (u) is a continuous func- 
tion of u. The Hamiltonian operators at the two times of 
measurement to and tf are selfadjoint operators. Accord- 
ing to the theorem of Stone [20] , each of the exponential 
operators exp {— iuH(to)} and exp {iuHu{tf)} forms a 
strongly continuous one-parameter group of unitary op- 
erators with parameter u. As the trace of a product of 
two strongly continuous operator valued functions of u 
with the density operator p(t ), which is a trace class op- 
erator and independent of u, the characteristic function 
(3) is a continuous function of u. 

Proof of property (ii): G tf .t ( u ) is a positive definite 
function of u. Using the cyclic invariance of the trace and 
the fact that H(to) and p(to) commute with each other, 
we can rewrite the left hand side of the inequality (7) as 



J2Gt f ,t (ui - uj)z*z 3 = ^-R e *(««-«i)**(t/) e'^ u '- u ^ H{t ^p{to)z*z 3 = TrA+Ap(t ) > , (Al) 



where 



E 



z n e —iv.iH H (tf) p iuiH(t a ) 



(A2) 



is a bounded operator and A + its adjoint. The last in- 
equality in (Al) immediately follows with the positivity 
of A + A and of the density matrix p(t ). 

Proof of property (iii): G f t,,t o (0) = 1. For u = the ex- 
ponential operators exp {— iuH(to)} and exp {iuHn (tf)} 
become unity. The trace over the density matrix p(to) re- 
duces by means of eqs. (5), (6) to the trace of the initial 
density matrix p(to), which is one. 



APPENDIX B: THE MATRIX ELEMENT 

(n| exp{iuH H (tf)}\n) 

The total time rate of change of the Hamiltonian Hn(t) 
coincides with its partial derivative with respect to the 
time which for the driven oscillator becomes, cf. eq. (17), 



dH H {t) 
dt 



r(t)a H (t)-f(t)a+(t), 



(Bl) 



r 



where an{t) and a#(t) denote annihilation and creation 
operators, respectively, in the Heisenberg picture, which 
are given by 



(I II' I i = e " "a 



a+(t) = e luJt a+ + % - [ dse^'^ f*(s 
n Jo 



(B2) 



This yields for H H (t f ) 

H H (t f ) = hiua+a + B*(t f )a + B{t f )a+ + C{t), (B3) 
where 

B(tf)= [ f dsf(s)e^ s 
Jo 

C{tf) = lf S dsj ds'\.f{s)f*{s')e lu{ - s - s '^ ( B4 ) 
n Jo Jo L 

-r( S )./V)e-^- s '>i 

The unitary operator 



V = e 



za^—z a 



(B5) 
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with 

transforms Hn{t) into 

VH H (t f )V+ = hwa+a + L(t f ), 

where 

\B(t f )\ 2 |/(i/)| 2 



L{tf) = C(t f ) 



(B6) 



(B7) 



(B8) 



Note that V induces a shift of the creation and annihila- 
tion operators 



VaV+ =a-z, Va+V+ = a+ - 



(B9) 



and further note that, when acting on the groundstate |0) 
with a|0) = 0, the operator V yields the coherent state 
\z), i.e. 



V\0) = \z). 
One finds with these properties 
1 



(BIO) 



(n|e 



iuH H (tf) 



n) 



,(z\(a-zr 



x e iuhuja + a+iuL(t f ) _ 



1 f) 2n 
1 c iuL(tA 



dx n dy n 



(z\e 



x(a — z)iuhLoa + ay(a + —z*) 



Z)\ X =y=0 

(Bll) 



Here we have introduced the auxiliary variables x and y 
which allow to represent the nth powers of shifted cre- 
ation and annihilation operators by derivatives of respec- 
tive order. The scalar function e ~ l i xz +v z ) can b e taken 
out of the scalar product and the remaining operator can 
be brought into normal order. It then becomes [23] 



e xa e »fi,«+a e9 «+ = ^ j exp ^ e iuhu> _ 1 j a+a 

(xa + ya++xy)]} , 



(B12) 



where under the normal ordering operator Af all creation 
operators stand left of the annihilation operators. The 
matrix element with respect to the coherent state \z) can 
be read off, yielding, 



(n\e 



iuH H (t f ) 



n) 



1 



e »£(t/)exp{(e' 



iHlo 



!) VA 2 } exp { (e^» - 1) (xz + yz*) + e^xy} \ x=v=0 



e »i(*/) exp |( e ^ _ y | z |2| __ ^ e iuhu, _ i) z + e ™^ y ] n exp{(e m?iw - 1) (\z\ 2 +yz*)} | w=0 



iu\f(t f )\ 2 /(hu;) 



exp{( e ^-l)|z| 2 }^ 



n\ \z 



2(n-k) 



iuhuk ( iuhu: 



\kj (n-k)\ 



1) 



2(n-fe) 



(B13) 



APPENDIX C: WORK DISTRIBUTION FOR A 
CANONICAL INITIAL STATE 

To determine the expression (36) for the work distri- 
bution q c (/3) we start from the general expression given 
in the first line of eq. (27). Interchanging the summation 
over the indices n and k we obtain 



oo Ik 



?m = e - |z| £ £(-iy 



m,k=0 1=0 



| z |2(fc+m) / 2k 

m\k\ ( I 



— r 



G) e -N 2 



fc 

~ j^L |_|2(fc+m) /oj. 

y V(-i)^ ( 2/c 

, / — n / — n \ 



-e-f> \k 



m,k=0 1=0 



lP -1 



( _i)- e -i-i 2 v Hz|2)m 



m=0 



mi 



x E 

fe=|m— r\ 



^-i 



(3) 



fc! 

(-w 2 r 



2fc 

fc + r — m 



2\z\ 2 /(e' 3 ~l) 



X L, 



m=0 

M 2 



ml 



11 



(±) e -|^| 2 coth(^/2) 



E 

m=0 



\z\ 2m ( 2\z 

i ^ |m — r| 



m! 



0-1 



(A) -|,| 2 coth(,3/2) (v^Pf ( 2 M 2 

1 ^ m! r ~ m U^-l 

^ m— 

oo I 

+ E - 



m=r-\-l 
~Ir—i 



Im—r 



ML 

e? - 1 



.e^-l 

e -| Z | 2 coth( / 3/2) e /3r/2 /j 



sinh(/3/2) 



(CI) 



In the first step ( = ) we performed the sum on n accord- 
ing to 



- - x K In 



Er 

n—k 



x\k 



1-x 



(C2) 



cf. Rcf. [24], 5.2.11.3. In the second step = the Kro- 
necker delta is used to perform the sum over k. The 

(3) 

third step = is based on the relation 

ii(^h e "'^ x> <C3) 

k—\l\ 

valid for integer I. Here I v {x) denotes the modified Bessel 
function of the first kind of order v. With I v {— x) = 



(—iyi v (x) where v is an integer, we come to the right 

(4) 

hand side of the equality = . In the next step the sum on 
m is rewritten. The term in the square brackets vanishes 
because I v {x) is an even function of the order v. The 
remaining sum can be performed by means of the identity 



00 fk /cy, \v/2 

fe=0 ' \ • / 



(C4) 



cf. [24] 5.8.3.1. This leads to the final result given in cq. 
(36). 



APPENDIX D: WORK DISTRIBUTION FOR A 
COHERENT INITIAL STATE EQ. (45) 

Starting from eq. (27) we may proceed in an analogous 
way as in the case of a canonical initial state, cf. the 
Appendix C. According to eq. (43) a Poissonian average 
over the binomial (^) has to be performed instead of the 
geometric average in the first step of eq. (CI). This 
yields 



kl 

n—k 



-\a\< 



|2fe 



k\ 



(Dl) 



Next the Kronnccker delta is used to perform the sum 
over / leaving one with two sums of which the inner one 
over k can be expressed in terms of a generalized hyper- 
geometric function, [21], to become 



E 

k—\rn — r\ 



(-M 2 ) 



2\k 



2k 



(-\az\ 2 )\ m ~ r \ ( 1 
■ \, ' - n IM2 1^2 |m-r| + -;|m-r| + l,2|m-r| + l;-4|az| 2 

(fc!)^ \k + r — mj (\m — r\\j z \ 2 



(D2) 



This immediately leads to the expression in eq. (45). 
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